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The latest results from the ATLAS and CMS experiments at the CERN Large Hadron Collider 
(LHC) unequivocally confirm the existence of a resonance, X, with mass near 125 GeV which 
could be the Higgs boson of the Standard Model. Measuring the properties (quantum numbers and 
couplings) of this resonance is of paramount importance. Initial analyses by the LHC collaborations 
disfavor specific alternative benchmark hypotheses, e.g. pure pseudoscalars or gravitons. However, 
this is just the first step in a long-term program of detailed measurements. We consider the most 
general set of operators in the decay channels X — > ZZ, WW, Z-y, 77 and derive the constraint 
implied by the measured rate. This allows us to provide a useful parametrization of the orthogonal 
independent Higgs coupling degrees of freedom as coordinates on a suitably defined sphere. 
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Introduction. The determination of the spin J and 
parity P of the putative Higgs boson [1] is of paramount 
importance. The observed decay to 77 [2] already ex- 
cludes the case of J = 1 by the Landau- Yang theo- 
rem [3] . Spin discrimination and coupling measurements 
are achieved by studying kinematic correlations of the 
final state objects, generally Icptons. An optimal vifay 
of incorporating all available information is through the 
Matrix Element Method (MEM) [1H7|. Initial results 
from ATLAS |8J and CMS HO] strongly disfavor 
J^ = 0~, as well as J^ — 2+ with graviton-like cou- 
plings. These results represent an important first foray 
into the general parameter space of generic boson cou- 
plings. 

In this paper we provide a theoretical framework for 
spin and coupling measurements in full generality, with- 
out theoretical prejudice towards specific benchmarks. 
For concreteness and simplicity, we consider the example 
of a spin-zero "X" resonance, decaying to four Icptons 
through two intermediate Z-bosons [TTI - fT^ . though our 
approach can be readily extended to other spins and final 
states. The complete measurement of the X couplings 
in full generality will provide insights into the nature of 
electroweak symmetry breaking and may offer the first 
glimpses of physics beyond the Standard Model (SM). 

Effective Theory. In what follows, we will consider 
a spin-zero state X. In general, X has no definite CP 
properties, and can be thought of as a linear combination 
of some CP-even state H and some CP-odd state A [15] : 



X = H cos a ~ A sin < 



(1) 



In the special case of a = 0, X is a pure J^ = 0"*" scalar, 
as predicted in the Standard Model, while for a — ±7r/2. 



X is a pure pseudoscalar with J^ — 0^ . 

The X couplings to two gauge bosons (e.g. ZZ) can 
be classified according to their symmetries in the CP- 
eigenstate basis {H,A). Three types of terms are possi- 
ble: 



Jlf'2 1 



(2) 



where P'"' is the Z-boson field strength tensor, Mz — 
91.1876 GeV is the Z-boson mass, v = 246 GeV is 
the electroweak scale, and /'■^•' and /^'*-' are (in gen- 
eral momentum-dependent) form-factors, which from the 
effective theory point of view should be thought of as infi- 
nite series expansions in terms of some new physics scale 
A. Let us discuss each one in turn. 

The form-factor flu, describes interactions of the CP 
even component H that necessarily violate SU{2) gauge 
invariance. Expanding in powers of A^^, one obtains 
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where gi are dimensionless coupling constants, and the 
derivative operator d {d) acts on the Z field to its left 
(right). Eq. ^ docs not include terms containing d^Z^, 
which vanish in the Lorenz gauge. 

The second class of CP-even couplings, described by 
the form-factor 
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may respect SU{2) gauge invariance, if iJ is a SM singlet. 
Note that terms of this sort could be generated by a 
[dfjdi, — g^iud''dpj term in Eq. (3|, which explains the 

absense of such a term there. Fmally, the couplings of 
the CP-odd component A are given by 
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and may be SU{2) gauge invariant as well, if A is a SM 
singlet. However, here we do not make any assumptions 
about the SU{2) quantum numbers of H and A; the 
terms Q and (|5| will also violate SU{2) gauge invari- 
ance if if or A (respectively) transform under SU{2). 

We note that the g^ and ge terms in ([3| are typically 
omitted in the literature. The contributions to the X — > 
ZZ amplitude resulting from the D + D and d^dp terms 
are proportional to Af|^ + Af|^ and M|- - M|^ - A/|^ , 
respectively, where Mx is the mass of the resonance and 
Mz^ and M^^ are the invariant masses of the two Z- 
bosons (in the usual convention where Mz^ > Af^^). 
Thus, if both Z-bosons were on-shell, the contribution 
from such terms would be essentially constant, and there- 
fore could be absorbed into a redefinition of gi. However, 
in the case of interest where Mx ~ 125 GeV, one or both 
of the Z's are off-shell and Mz^ and M^^ vary from event 
to event. Thus, strictly speaking, the g^ and ge terms 
cannot be absorbed into gi, though their effects (rela- 
tive to gi) are expected to be rather small, due to the 
additional A^ suppression. 

It appears that a general analysis of the couplings 
of the Higgs boson candidate would have to include 
at the very minimum the terms identified in eqs. pi 
pi), and perhaps even the higher dimensional operators 
~ A~^,A~^,... which were not explicitly listed. Such 
an analysis may indeed be desirable at some point in 
the future, when significantly more data will have been 
accumulated by the LHC experiments. However, there 
are strong theoretical and experimental motivations for 
making certain simplifying assumptions for analyses in 
the immediate future. 

First, all experimental evidence so far suggests that the 
new physics scale A is high compared to v. Second, con- 
sistency of the effective theory description requires that 
A be sufficiently removed from the relevant experimental 
energy scale. Finally, the higher dimensional nature of 
the 1/A couplings suggests a radiative origin and hence 
suppressions by loop factors. It is therefore reasonable 
to expect that the higher order terms 53, 175, 176, etc. in 
the expansions pfsl) are negligible relative to the corre- 
sponding leading order couplings gi, (72 and (74. 

At the same time, the relative size of the leading terms 
gi, g2 and 54 is a priori unknown. For example, the 52 
and ^4 terms are equally suppressed by A, and may both 
preserve gauge invariance, thus it is difficult to argue that 
one should be preferred over the other. Similarly, we do 



FIG. 1: The principle of geolocating: the change of variables 
( 9 1 inflates (a) the hyper-surface defined by the measurement 
(81 into (b) a perfect sphere. 



not know the extent to which H is involved in breaking 
SU{2) gauge invariance, and hence we cannot simply as- 
sume that the SU{2) breaking term gi dominates over 
52 , in spite of the additional A suppression in the latter 
[T^ . This is why we keep all three terms and consider 
the effective Lagrangian for the mass eigenstate X to be 



C = X 



Kl- 



^z^z^ 



K2 

2v 



Fp^F^- 



K3 
2v 



Fp.F^' 



, (6) 



where 



Ki=giCOsa, K,2 = g2-j^ cosa, n^ = gn^—sma (7) 
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are the effective couplings which need to be measured by 
experiment. We note that this procedure means neglect- 
ing (for the reasons stated above) the (73 operator from 
Refs. [F, "5" and the 95 and gg terms listed in (|3| . 

Counting degrees of freedom. From an effective 
theory point of view, the coefficients k; in the Lagrangian 
([6]) are real numbers. Thus the problem of measuring the 
Higgs couplings is reduced to determining the three di- 
mensionless degrees of freedom kj. However, one can 
go one step further and use the measured signal rate 
{X production cross-section times branching fraction for 
X — > ZZ) to eliminate one more degree of freedom, 
namely the overall k scale. In order to do this, we will 
assume that the X can be treated as a narrow resonance. 
We then note that the partial width for X -^ ZZ is a 
quadratic function of the k's: 



T{X ^ ZZ) = TsM Y, 
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where we have factored out the partial width T sm pre- 
dicted in the SM in order to define dimensionless coef- 
ficients 7y listed in Table IT| The measured total rate 
then provides the overall normalization and constrains 
the Ki couplings to lie on the closed hyper-surface shown 
in Fig. fly a). The idea now is to change variables and 
parameterize the couplings Ki in terms of the two coor- 
dinates on this hyper-surface. This re-parametrization 
is useful and meaningful, since the (normalized) angu- 
lar and invariant mass distributions used to measure the 



TABLE I: Numerical values for the coefficients defined in ([8| 
for different diboson final states. 



K, Map 



J Map 



Process 


7ii 


722 


733 


7l2 


X^ ZZ (DF) 
X^ ZZ (SF) 

X->77 
X^ WW 


1 
1 

1 


0.090 
0.081 

1 
0.202 


0.038 
0.032 

1 
0.084 


-0.250 
-0.243 


-0.379 


after cuts 


X ^ ZZ (DF) 


1 


0.101 


0.037 


-0.277 



Higgs spin and CP properties are insensitive to the over- 
all scale of the couphngs k^. OperationaUy, we propose 
to do this by changing variables as 



7 . Oij Kj , 
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where O21 =031 = O32 = and 

Ou = lu/Vl^, {i = 1, 2, 3) 
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(z = 2,3), 



In terms of the new variables Xi, the constraint implied 
by (Is]) is the surface of a sphere, as seen in Fig. [ijb). 

Note that in the case of real couplings, 713 and 723 
vanish identically (and hence are not listed in (Is])). This 
can be understood in terms of parity; interference terms 
between amplitudes describing the decay of a parity-even 
scalar and a parity-odd scalar are odd under parity, hence 
vanish under integration provided the cuts respect par- 
ity (as is nearly always the case). The situation in the 
presence of cuts is less clear if the couplings of the three 
operators in ([6]) are complex valued; however, any 713 or 
723 term generated in this case is expected to be very 
small. 

Geolocating the Higgs. A sphere may be pa- 
rameterized using latitude (p and longitude A. Making 
an obvious analogy, we can think of the tree-level SM 
Higgs couplings (ki,K2,K3) = (1,0,0) as being repre- 
sented by the point (0, A) — (0, 0) in the Gulf of Guinea 
(see Fig. [2|. The case of a pure pseudoscalar with 
(ki, K2, K3) = (0, 0, 1) corresponds to the North pole with 
{(j>, A) = (90, 0). Those are the two main benchmark sce- 
narios considered so far by the LHC experiments. The 
Higgs sphere from Fig. WCh) now opens up the full range 
of possibilities. For example, one may consider a cer- 
tain amount of mixing as in (fTl) between the SM Higgs 
and a pseudoscalar, placing us on the Greenwich merid- 
ian. Alternatively, one may allow non-trivial values for 
the two CP-even couplings ki and K2, thus spanning the 
equator. Finally, one could also envision the most gen- 
eral case with non-trivial values for all three couplings. 
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FIG. 2: The effective couplings Ki and the signal efficiency 
after cuts as a function of latitude (0) and longitude (A) for 
the different flavor 4 lepton channel X — ?> ZZ — >■ 2e2p. 



e.g. {ki,K2,K3) = (-0.945804,-3.88525,2.44522), cor- 
responding to (0,A) = (29.64945,-82.3486), which hap- 
pens to be a location in the south end zone of the Swamp. 
We note that CMS has recently undertaken [TD] an 
exploration of the sphere along the Greenwich meridian 
(A = 0) by measuring the parameter fa3 = 1^3^/(1^3^-1- 



|j4ip), parametrizing the relative contributions |^ip and 

nch 
(0,0,1) 



IA3P to the total cross-section of the two benchmark 



models {ni, K2, 1^3) = (liO,0) and (ki,K2,K3) 
respectively. 

Selection cut bias. By definition, the measurement 
of the Higgs signal rate (and, by association, of the par- 
tial width ^) is done using only events which pass se- 
lection cuts. Therefore, it needs to be corrected for the 
efficiency. Unfortunately, as demonstrated in the bottom 
right panel in Fig. [2| the signal efficiency is noticeably 
model dependent, i.e. it is quite sensitive to the actual 
values of the individual Ki parameters, and may vary 
from as low as 35% to as high as 54%. One source of 
this variation is the cut on M^^ (see e.g. [H]). The ef- 
fect is quantified in Fig. [31 which shows unit-normalized 
Mz2 distributions for several choices of Ki, including an 
example with a high efficiency {{ki, K2, K3) — (0,0,1), 
black circles) and a suitably chosen example with a low 
efficiency ((ki,K2,K3) = (1.8,5,0.08), red crosses). The 
Mz2 distribution in the latter case is softer, and many 
events fall below the minimum accepted Af^^, value of 
12 GeV. We conclude that it would be virtually impossi- 
ble to correct for the efficiency without knowledge of the 
couplings Ki, which we are trying to measure in the first 
place, thus falling into a vicious circle. 

Here we advocate an alternative approach. It is suf- 
ficient to realize that the coefficients jij, while also af- 
fected by the efficiency, do not depend on Ki. Therefore, 
in defining the Higgs sphere ^ one could simply use the 
corresponding values of jij after cuts, which need to be 
calculated once and for all. As shown in Table |l) the 
changes are subtle, yet noticeable. Note that even with- 
out any cuts, the same flavor (SF) and different flavor 
(DF) coefficients are different, due to the interference ef- 
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FIG. 3: Unit-normalized Mz2 distributions for four different 
tlieory bendimarks (Ki,fi:2,K3). 



(K„K^,K,)=(1,0,0) 




. (K„K„K3)=(0,1 0) 



^•1-^ — fr^ — % 



J 



longitude 



J 



-1 

longitude 



FIG. 4: Distribution of tlie best fit values for the Higgs 
geolocation in 1000 pseudoexperiments with 300 X — >■ ZZ -> 
2/i2e signal events each, for four different benchmark points. 



fects present in the SF case (see, e.g. [6]). 

Interpretation. We emphasize that by studying the 
angular and invariant mass distributions of the final state 
leptons in the X -^ ZZ ^ At channel it is experimen- 
tally possible to determine the exact geolocation of the 
Higgs boson candidate, without any simplifying theoret- 
ical assumptions. A proof of principle is shown in Fig. [4] 
We consider 4 benchmark scenarios for {ki, k,2t nz)'- (A) 
(1,0,0), (o) (0,1,0), (*) (0,0,1), and (+) (1,1,1). We 
show results from 1000 pseudoexperiments with 300 sig- 
nal events each (after cuts). Fig. H shows that in each 
case, the maximum likelihood fit indeed selects the cor- 
rect geolocation (marked with an open symbol) of the 
Higgs candidate. 

(In)deterniination of the mixing angle a. We 
note that geolocating the Higgs candidate by itself is not 
sufficient to determine the value of the mixing angle a, 
and we will still not know the relative composition of X 
in terms of H and A. As seen in eq. (l7|, the measurement 
of all three effective couplings ki, K2 and K3 only places 
3 constraints on the 4 input parameters a, gi, (72/A and 
(74/A, so one degree of freedom (which can be chosen to 
be the mixing angle a) will always be left undetermined. 

Music of the Spheres. Until now, we have been 



treating the coefficients Ki of the operators in the effec- 
tive Lagrangian (|6| to be real. We term this "Scenario 0" . 
Of course, this must be the case at tree level, and even be- 
yond, as long as those operators are generated by loop di- 
agrams involving heavy particles, such as the top, which 
may be consistently integrated out of the full theory to 
form an effective theory. However, in accordance with 
the optical theorem, loops with light particles (such as 
&-quarks, which cannot be integrated out of the theory) 
will contribute imaginary parts to those couplings. In 
this case, each of the ki is complex and there are five 
physical degrees of freedom (an overall phase in the Hi 
may be removed). Nevertheless, the width constraint 
T{X — >• ZZ) = ^SAi'^i jJij^iKj can still be enforced. 
As before, the diagonalization and rescaling (|9| renders 
this a sphere in C'^ (rather than an ellipsoid, as is the 
case in general). 

Even in this general complexified case, there are three 
interesting and simple scenarios where the number of de- 
grees of freedom is reduced to three. Namely, if one of Ki 
is too small to have observable effects, then the remain- 
ing two Ki, taken to be arbitrary complex numbers, also 
comprise two observable degrees of freedom (as an over- 
all phase may be removed). There are, of course three 
such scenarios, "Scenario 1" where ki is negligible, "Sce- 
nario 2" where K2 is negligible, and "Scenario 3" where 
K3 is negligible. 

All four scenarios have physical motivation. Scenario 
holds when the important couplings occur at tree-level 
or come largely from top-loops. Scenario 1 describes the 
most general couplings to a massless particle such as the 
photon or the gluon, though it could be obtained in the 
ZZ or WW cases if the X has nothing to do with EWSB. 
Scenario 2 describes, e.g., a state where X is primarily 
A, but the H with which it mixes is an SM-like Higgs bo- 
son. Scenario 3 describes AT as a general CP-even scalar. 
The scenarios are not totally exclusive; for instance, the 
ZZ case in the Standard Model is well approximated by 
points in Scenarios 0, 2, or 3. 

Other final state channels. While we have so far 
only focused on the ZZ — > 4£ final state, the preceding 
discussion can be readily applied to the other di-boson 
final states, W^W' , Z7, and 77. The interactions rel- 
evant for those channels involve corresponding new cou- 
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and t^P ■ Assuming SU(2) gauge sym- 



metry, these new couplings are related, e.g. 
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with obvious definitions of /t's. As before, the measured 
rate in each of those channels defines a sphere similar 
to eq. (|8|, where the corresponding values for the coef- 
ficients 7ij are given in Table ll| It is interesting to note 



that these spheres are not overlapping (even for the DF 
and SF case in the same X — >■ ZZ channel), and in princi- 
ple couplings can be measured from intersecting spheres, 
although we expect that this method is of only academic 
interest, since the corresponding precision would be very 
poor. 

Experimental implementation. While the exam- 
ple analyses performed here are simplified, the additional 
steps needed to perform an actual analysis are relatively 
straightforward. One needs to 

1. Calculate the parameters ^ij from Table IT] using 
the actual detector geometry and effciencies and 
the specific cuts to be used in the analysis. 

2. Use the 7ij to define a sphere in coupling space, 
either in the scenario where all couplings are real 
(Eq. (|8|) or in one of the other three scenarios. 

3. Determine the point on the sphere in coupling space 
that maximizes the likelihood associated with the 
observed events. The likelihood employed should 
include background information f5H7tll4j. 

One potential experimental concern stems from the 
difhculty of generating events to perform pseudo- 
experiments for sufficiently many points in the contin- 
uous parameter space. In experiments, the significance 
of a measurement is often determined by comparing the 
observed value of some test statistic to a distribution of 
that test statistic obtained from pseudo-experiments. In 
principle, such an approach is not necessary when the 
test statistic used is the likelihood, since the likelihood 
has a well-defined statistical meaning. However, one may 
still wish to perform pseudo-experiments in this case, for 
example to take into account reducible backgrounds not 
included in the likelihood expression used. 

We point out that this concern is unjustified; one can 
perform pseudo-experiments without generating events 
for every point in parameter space. Our recommended 
procedure is to generate events for one particular point 
in parameter space, which we label by ((/)o,Ao). Then 
one can perform pseudo-experiments for a different point 
in parameter space, {(ptestT^test), using the same events, 
but calculating the likelihood not from the usual formula 



ln£(0,A) = ^lnP(p„</),A), 



(13) 



but from 
ln/:(</>,A) 



"(Pij Vtest, \est) 
P(Pj,0O,Ao) 



lnP(p„0,A), (14) 



where Pi is the kinematic information in the i-th event. 
In other words, we suggest recycling the events in the 
pseudo-experiments by reweighting, which is relatively 
trivial by using JHUGen [F or MEKD [6 , rather than gener- 
ating a new event sample for each parameter space point. 



Summary. We have proposed several related parame- 
terizations of the couplings of the 125 GeV boson discov- 
ered at the LHC. These parameterizations allow the LHC 
experiments to go beyond comparing benchmark points, 
etc. in a relatively simple, yet very general, theoretically 
motivated fashion. The low dimensionality of the pa- 
rameter space is helpful experimentally as well as in the 
visualization and interpretation of results. We look for- 
ward to the implementation of this method by the LHC 
collaborations and to the continued exploration of the 
couplings of the putative Higgs boson. 
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